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ABSTRACT 


A horizontal cylinder, located near the floor el 7am o 
dimensional wave channel was subjected to a train of non- 
linear gravity waves. The horizontal and vertical forces 
were measured and presented in dimensionless form. Experi- 
mental values of horizontal and vertical force coefficients 
are presented as functions of dimensionless wave height and 
dimensionless wave period. The dimensionless force coeffi- 
@#ents predicted by a modified Morrison equation and a 
Froude-Krilov force are compared to experimental data. 

Fluid particle velocity and acceleration values were calcu- 
lated erom Stokes fifth-order gravity wave theory. Experi- 
mental dimensionless wave periods from 10 toeecO0 were 
investigated. 

The horizontal force coefficients were found to vary 
linearly with wave height for dimensionless period values 
from 60 to 120. The vertical force coefficients were 
found to be inertia dominated at low dimensionless wave 
heights and dominated by a lift force at higher wave heights. 
The theory predicted experimental values of force coefficients 


with good accuracy, especially at greater water depths. 
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ieee CRODUCT ION 


Exploration for and use of ocean resources has received 
increasing attention in the scientific community in recent 
years. A considerable effort has resulted in the development 
of technology in areas such as off-shore oil exploration, 
recovery of petroleum and deployment of underwater habitats. 
Oil production often involved the use of very large submerged 
oil storage tanks and associated piping systems. In some 
locations, oil is transferred long distances in submerged 
piping systems, often leading through the surf zone to 
facilities ashore. Also, waste disposal systems usually 
involve deployment of large diameter outfall pipes which are 
laid through the surf zone. Proposals for large scale mining 
and food gathering activities indicate the trend toward 
design and construction of large submerged structures. 

This activity has generated considerable interest in the 
interaction of gravity waves with submerged objects and 
particularly in the forces induced. Most previous investi- 
gations [1,2,3] in this area center on the use of a so-called 
Morrison Equation [3], which involves both a drag and iner- 
tial component of force. Unfortunately, the drag and added 
mass coefficients are found to be neither constant nor 
simply shape - or Reynolds number - dependent in the unsteady 
flow that occurs in the interaction of waves with fixed 
bodies. In addition, where the amplitude of wave motion is 


large compared to the body dimensions, flow separation occurs, 





and the coefficients vary as the flow in a given direction 
develops. Finally, for small depth to characteristic dimen- 
sion ratios, the proximity of the free surface may influence 
the value of the coefficients, even if no flow separation 
Seeurs. 

In spite of these difficulties, the Morrison Equation 
will, with certain restrictions, provide a complete descrip- 
Bion of the horizontal forces acting on a large horizontal 
cylinder in contact with a plane boundary. If the cylinder 
is large in comparison to the amplitude of the fluid motion, 
separation does not occur, and the resulting wave-induced 
oscillatory flow may be considered an unséparated potential 
flow about the cylinder. This has been verified by Kuelegan 
and Carpenter [2] and by Sarpkaya and Garrison [4]. For 
such cases, the drag component of force is negligible, while 
the added mass coefficient is a constant. 

Application of the Morrison Equation to determine the 
vertical components of force for a bottom-mounted cylinder 
yields a vertical force component of zero, an obviously 
imvalid result. It is recognized that the vertical force 
should be composed of two components; one associated with 
uae inertia.effectseand a second which accounts for the lift 
force caused by the increased velocity and hence decreased 
the pressure over the top of the cylinder. 

The inertial component of vertical force can be approxi- 
mated by what is commonly known as a Froude-Krilov force. 


This is the force caused by the pressure distribution around 


10 





the surface of the cylinder existing in the wave if the 
Cyiander was not present. This is an approximagdonmm iach 
does not account for the presence of the body and is 
therefore considered to be an underestimate of the force. 

The lift component of vertical force can be approximated 
by use of the unseparated, potential flow model used in the 
analysis of the horizontal forces, i.e., that of a uniform 
flow past a cylinder in contact with a rigid boundary. The 
lift coefficient for this case has been determined by Dalton 
and Helfenstein [5]. 

A study carried out by Johnson [6] for horizontal forces 
acting on a bottom-mounted, horizontal cylinder in long waves, 
eesumed that at long wave lengths, the horizontal forces 
depended on wave height and water depth only. This assumption 
is, however, of doubtful validity in the range of Johnson's 
west. 

Shiller [7] measured wave forces on a submerged horizontal 
cylinder due to small amplitude waves. Over the range of 
wave heights considered, the magnitude of the horizontal 
ieee was found to vary linearly with wave height. The 
magnitude of the vertical force was found to increase in 
proportion to the wave height squared. 

Perkinson [8] extended Shiller's investigation to include 
larger wave heights and periods. His study indicates that 
horizontal forces increased linearly with wave heights and 
becomes independent of wave period for large wave periods. 


The vertical forces were found to contain two regimes: the 


diet 





lower wave lengths where lift force predominates; and the 
Paeher wave lengths where inertial force is deminarre 

In order to evaluate any forces using the above analyses, 
it is necessary to first evaluate fluid velocities, accelera- 
tions and pressure by use of some wave theory. Fairly simple 
expressions may be derived for the forces by use of Airy wave 
theory. However, for finite amplitude waves in finite depth 
Wecer, a non-linear wave theory is considered to be necessary. 
stokes fifth-order wave theory has been chosen for this 
siudy . 

Stokes originally developed a second-order theory [9] for 
the case of a non-linear wave in water of finite depth. This 
method has been extended by Borgman and Chappelear [10] to 
third order. Skjelbreia and Hendrickson [11] have extended 
the solution to third and fifth order. Bretschneider [12] 
Mes presented a method for extension to any order. 

The purpose of this study is to compare the analytical 
results of Stoke's fifth-order wave theory with experimental 
values at longer wave periods. The transition section of the 
wave channel used by Perkinson was extended approximately 
forty feet ae provide a slope of 1:20. This provided the 
possibility of extending Perkinson's work to include longer 
meeyve periods. This study, then, is an extension of the 


work of Shiller [7] and Perkinson [8]. 
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II. THEORETICAL ANALYSIS 


A. PROBLEM DEFINITION 

The problem under consideration is illustrated in Figure 
1. A train of gravity waves is considered to progress in 
the positive x-direction in water of depth h, with rate of 
propagation C. The fluid particle velocity is expressed in 
terms of the horizontal and vertical components, u and o. It 
is of primary interest to determine the horizontal and 
vertical components of wave force acting on the horizontal 


cylinder which is in contact with the rigid bottom. 


B. DIMENSIONAL ANALYSIS 

An exact analytical solution to this problem is quite 
formidable and, therefore, an approximate solution is con- 
sidered in this work. However, first it is instructive to 
carry out a dimensional analysis of all pertinent parameters 
movolved. 

In general, the maximum of the wave force per unit 
length acting on a submerged cylinder due to wave motion is 


known a priori to depend on the following parameters: 


ih ae (Se Pt in) OnE oly Gin) 
where 

F = wave force 

h = water depth 

H = wave height 

T = wave period 

L = wave length 


1s ; 





T= Wave Period 





Fig. |. DEFINITION OF GEOMETRY 
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= cylinder radius 

= fluid density 
gravitational constant 
= fluid viscosity 


2 Fr TM UD Pf 
tt 


= cylinder length 


However, a relationship exists between the parameters 
associated with the incident wave (i.e. h,H,T and L). | 
Consequently, only three of these parameters (e.g. h,H and 
T) are needed to completely specify the incident wave. A 


dimensional analysis of the parameters in equation (1) yields 


the following groups: 


F /ogaé 2 = f.(gT/n, h/as gy ane 7 & V pha‘) (2> 


max 


The last term on the right side of equation (2) repre- 
sents the ratio of Froude number to Reynolds number, indi- 
@eacving the ratio of viscous to inertial forces. It is 
believed that, if this number is small, it may be neglected 


from further consideration, and equation (2) may be rewritten 


oe ie 2 me oor 
ea = f2(et /n, H, h) (3) 
where F = F /0 ack 
max max’ ?5 
H = H/2a 
h = h/a 


The approximate analytical approach of this study seeks: 
1. To develop expressions for the fluid particle motion 
in the gravity wave as functions of the incident wave 


parameters, and 
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2). To develop expressions for the forces on the sub- 
merged cylinder as functions of the fluid particle 
motion. 

The assumptions upon which this analysis is based are: 


1). Water will be considered inviscid and incompressible, 


2). The cylinder radius will be considered small in 
comparison to the dimensions (h and L) of the 


incident wave. 


C. WAVE THEORY 
The motion of fluid particles in a gravity wave based on 
assumptions 1) and 2) above is specified by the following 
boundary value problem. 
1). Governing Equation 
Assumption 1) implies that Laplace's Equation 


governs the fluid motion. That is: 
veo = 0 (4a) 


where $¢ is the velocity potential and the two velocity 


components are given in terms of @¢ as 


an: 
» ve 7 (4b) 


- 9¢ 
: oo 

2). Boundary Conditions 
The boundary condition on the bottom is that of a 


non-porous wall, i.e., 


1200 


ay = Sie y = —h (5) 
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Two boundary conditions are encountered at the water 
surface, one dynamic and one kinematic. 

The dynamic boundary condition is obtained from 
Bernoulli's Equation by setting the pressure equal to zero. 


imac is: 


ue + ve + 260= 2e(0 - y) at y=n (6) 


where Ii is the total energy head, a constant. 

Assuming the wave travels without change in form, it is 
then possible to choose a reference system moving in the 
positive x direction with wave celerity C. This makes the 
fluid motion steady with respect to the moving reference 
system. The Bernoulli Equation in such a reference system 


es: 


(u-C)© +v%= 2e(l-y) atyen (7) 


The kinetic boundary condition imposes the condition that 
no fluid be transported across the free surface. Or, in 
other words, the free surface must be a streamline. This may 


be @xpressed mathematically in the form: 


2 


ui x 
i WSC we me (8) 





Equations (4a), (5), (7), and (8) completely specify the 
boundary value problem. 

There are, historically, several approximate solutions 
to this problem. The solution to be considered further in 
this analysis is that commonly known as Stokes fifth-order 


wave theory. 
1 





Pee STOKES FIFTH ORDER SOLUTION 

Stokes proposed [9] a power series type of solution to 
the above problem and proved the validity of the second order 
solution based on this power series. Skjelbreia and 
Hendrickson [9] have extended this theory to fifth-order, 
yielding a solution to the boundary value problem of the 


form : 


cosh(28S)sin(26) 


2 : 
d = gi A, cosh(BS)sing + A, 


+ A,cosh(3BS)sin(36) 4 A,cosh(4gS)sin(46) 


+ A-cosh(58S)sin(56)} (9a) 


where Ay through Ag 


This equation may be written in dimensionless form as: 


are as given in Appendix A. 


cosh(28S)sin(26)+...} 


Ot 
| 


= C{A,cosh(8S)siné + A 


; (9b) 
where: $ = £6/\/ gh 
C = C/\/ gh 


Differentiation of equation (9a) with respect to time 
yields 


22 c“{A cosh(8S)cosé@ + 2A,cosh(28S)cos (286) 


dt 2 


+ 3A,cosh(38S)cos (38) s 4A,cosh(4BS)cos (46) 


+ 5A,cosh(58S)cos(58) } (10a) 


ab 





or, in dimensionless form 


~~ 


9@ = é*{ A, cosh(BS)cos(é) een 


aE cosh(28S)cos(26) 


@ 


4 
aa) 


(Oise) 
where 22 = oe /\Jen 


The series form of the wave profile is assumed to be: 


y = z{Acos (0) + B,cos(20) + B,cos(36) 


e 3 


+ B,cos(46) + B.cos (56) } (ics) 


where Bo Throusa Be are as given in Appendix A and iX is an 
Seoterary constant. 


Equation (lla) can be expressed in dimensionless form 


as 
y = a cos(é) + B,cos (26) + ...} (11b) 
qd 
where y = y/a 
qd = onmh/L 


The series form for II is assumed to be: 


+ ANC) I (12) 


- i74¢ 
N= FA°C, 


where C. and Cy are as given in Appendix A. 


The wave celerity expressed in series form is 


ae 2 2] 
6: a C= + °C, + AAC. } : (13a) 
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or, in dimensionless form 


2 2 
cr tl + A°C, + MHC} mca cies) 


ona y= 


where Cy and Co are as given in Appendix A. 


Using equations (4b) and (9b) we find the horizontal 


particle velocity to be: 


u = (32) = C{A, cosh(8S)cos(8) + 2A.cosh(2B8S)cos(26) 


2 


+ 3A,cosh(3BS )cos (36) “+ 4A, cosh(4BS)cos (46) 


+ 5A-cosh(58S)cos(56)} (14) 


> 


where u = u/\/ gh 


Also, from equations (4b) and (9a): 


v= (ee) = C{A,sinh(8S)cos(0) + 2A.sinh(28S)cos(26) 


2 


+ 3A.sinh (38S)cos(36) + 4A, sinh(4BS)cos (46) 


3 
+ 5A-sinh(58S)cos(50)} (15) 


where v= v/\Ven 
Differentiation of equations (14) and (15) with respect 


to time yields the particle accelerations 


Soe Oe dc*{a,sinh(BS)sin(8) + 2A 


x 5 sinh(28S)sin(26) 


2 


+ 3A.sinh(38S)sin(30) + 4A, sinh(48S)sin(46) 


3 


+ 5A-sinn(58S)sin(56) GES 


5 


where 2, = ay/ eB 
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and 


ay = d¢*{Aa,sinh(BS)cos(8) + 2A,sinh(2BS)cos(26) 
: | (17%) 
where a. = a /& 


The potential function and it's derivatives, as well as 
the wave profile are now known in terms of the unknown 
quantities A and d. Specifying equations for these parameters 
are found by developing equations for the peak-to-trough 
wave height and the wave period. That is, the dimensionless 


wave height is defined as: 


Hew Sagan ogee (18) 
where Hd = H/h | 
Substitution of equation (13b) into equation (18) 
yields: 
id =820) tea + 1? (B [Bo (19) 
5 Si 33 55 
where B33; Boos and Bes are as given in Appendix A. 


Now, the wave period and the wave celerity are related 


by the relationship: 
¢ ="b7E (20) 


So, equation (15a) becomes 


2 
C 

(L/T)* = g—“f1 + ig + r tc) (21) 
; | 


eal 





Or 5 


ww 


d 
TS ihe 


ae + 2c. + 27 


es 
h/gT = 1 





Co} (22) 


The parameter eT /h is referred to herein as the period 
Deamameter. 

Equations (19) and (22), if solved simultaneously, yield 
values for A and d in terms of Hd and eT“ /h, quantities which 
are incident wave parameters. Thus the potential function, 
particle velocities and accelerations are shown to be functions 
of the dimensionless parameters obtained by dimensional 
analysis and given by equation (3). (Note that Ha = 2H/h) 

ihe next step in the analysis is to express the forces 
eetane on the cylinder as functions of the potential Dune tater , 


particle velocities and accelerations. 


E. FORCES ON THE CYLINDER 
Mime forces acting om the cylinder are described in terms 
eseeeemeir horizontal and vertical components. 
1). Horizontal Component 
The horadzontal component of force acting on the 
cylinder is expressed in the form of the so-called Morrison 


Equation as. 


2 egies?) + (1.0 + C na“ La (23) 
x 2 . M On 
where Cy = drag coefficient 


Cu = added mass coefficient 


ae 





For cases where fluid particle motion is smaiivin 
eomparison with the cylinder diameter, flow Separactemeeo-c 
Mee occur and the contribution of drag to the total force 
may be disregarded. In this case, equation (23) may be 


rewritten 


My 
Ul 


CEO Cy ay (24) 


F/pga°n. 


where F 

The added mass coefficient for a circular cylinder in 
contact with a rigid wall was given in closed form by 
Garrison [13] as Cy = 2-29. 

2). Vertical Force Component 

' As mentioned previously, application of the Morrison 
BPomacion to the vertical compoment of force yields a result 
which is clearly erroneous. The vertical force may be 
@emcraaered to arise from two sources: an inertial force 
mee a lift force. 
a). Lift Force 


The: UehCeLorce May oe expressed ais - 


ny 


1/2(2a%p)C,ulul (25) 


where c life coefficient. 


L 


or, in dimensionless form: 


FS hc, ul ul | (26) 
where F =F Joga’ 
L L ° 
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Dalton and Helfstein indicate [7] that C,=4.49 
for unseparated flow past a circular cylinder in contact 
Mth a rigid wall. 

b). Inertial Force 

The pressure distribution around the surface of 
imme cylinder extsting if the cylinder is not present results 
in a net vertical force (known as a Froude-Krilov force). 

Consider the resultant differential force dF. 
acting on a differential area of the cylinder dA as shown 
in Figure 2. The above equation, written in dimensionless 


fmemm, is: 
Bo San P(a,y)sinydy (30) 
The pressure is obtained from the Bernoulli Equation: 

P= -p{5(ue + ve) + So) _ (3) 


or, in dimensionless form: 


i eee a¢ 
P = -h[5(u + vo) + GS) | (32) 
where P = P/pga 
Therefore 
~ a ee SO ry 
F, = -h s[3(ae + vo) + ($8) Jsinya¥ (33) 
0 


And the total vertical force may be written: 


2 Il = 
B h{C,ulu| : [=(u + vi) + a 
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Fig. 2 


PER TIAL PORCE GEOMETRY 
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F. DISCUSSION OF COMPUTER PROGRAM 

The equations developed above were programmed for a 
eeeital computer as given in Appendix D. The basic inpugs 
to the program are H, h and et /h. The program yields 
Memvticle velocity, acceleration and computes the horizontal 
and vertical forces on the cylinder at various wave angles. 

are WlAnbo) SubRPetitame cadculates the coefficients for 
the Stokes fifth order equations. Simultaneous solution 
to equations (19) and (22) are obtained by assuming an initial 
value of X and d based on linear wave theory and iterating 
to the actual value using the Newton-Raphson technique. 

The main program generates tables of values and in-line 
graphs of the wave profile, particle velocity vs. depth, 


vertical and horizontal force vs. phase angle. 
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Iii. DESCRIPTION OF APPAR At Sa ue 
EXPERIMENTAL PROCEDURE 


A. WAVE CHANNEL 

Wave forces on the cylinder were measured experimentally 
in a wave channel. The basic wave channel, as shown in 
Figure 3, was fifteen inches wide and consisted of two sec- 
tions: a transition section and a shallow water section. 

The overall channel length was 87 feet. 

Three-quarter inch exterior plywood sheets were used for 
Memes channel sides. Vertical rigidity was provided by bracing 
the plywood sheets with 2x4's at various intervals. The 
vertical members were placed four feet on centers in the 
Shallow water section and from four feet to eighteen inches 
im the transition section. Additional rigidity was provided 
in the transition section by truss arrangements, as shown 
in Figure 4. Two 2x4's were attached to the top of the 
Vertical studs and extended the length of the channel. 

A double floor constructed of two plywood sheets with 
two half-inch separators was used. The sides of the channel 
were bolted together, through the spacing of the double 
Meeecon, with three-eights inch threaded rod. 

All wooden sections exposed to water were water-proofed 
and painted with three coats of Morewear Vitri-Glaz 1320A. 
Dow Corning Sealant 780 was applied to all seams and joints. 

A paddle type wave generator was used to generate waves 


in the channel. An aluminum plate was hinged at the tank 
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Fig.3 WAVE CHANNEL 
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Fig.4 


TRANSITION SECTION TRUSS ARRG T. 
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fioor and attached by a pin connection at the top to a drivie. 
rod. Wave motion behind the paddie was damped by use of 
baffle plates. A two horse power variable speed drive, with 
an output speed range of twenty to one hundred eighty revolu- 
tions per minute, was mounted on top of the wave channel and 
fitted with a six and three quarter inch radius face plate. 
Ways were attached to the plate to adjust the driving rod 
eccentricity between zero and eight inches. 

The transition section contained a ramp with a slope of 
1:20 which shoaled the generated wave with a change in depth 
from four to two feet. The overall length of the transition 
section was 42 feet. 

The shallow water section had a depth of two feet and 
an overall length of 42 feet. The test module was located 
twenty four feet from the transition, allowing the waves to 
reach a fully developed state following the transition before 
reaching the test section. 

A variable slope dissipater beach was located at the end 
of the test section in order to dissipate the wave motion. 
The beach was constructed of metal shavings held between two 
pieces of perforated stainless steel sheet metal, which were 
separated by two inch wooden spacers. A solid sheet of one- 
eighth inch aluminum was attached to the bottom of the beach 
by three and five eighths inch separators. One by one inch 
aluminum angles were attached to the exposed surface of the 
beach, parallel to the wave fronts. Maximum slope of the 


beach was 1:7. 
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B. DESTeMODULE 

The test module consisted ef the circullarecsawaer with 
its. own floor and walls, constructed as an integral system 
to provide easy removal of the model from the tank. The wave 
channel sides and floor were recessed to provide a smooth 
Eratiort1on irom channeleto test module. The channelewalls 
were made of plexiglass in the test module region to allow 
maximum visibility of the model. The test module is shown 
in Figure on 

A four inch diameter plexiglass cylinder was supported 
between the walls of the test module by use of cantilever 
beams and adjusted to approximately one-sixteenth inch above 
tne floor. A thin flexible plastic barrier was installed in 
miemeeap tO prevent any water motion under theseylinder due 
to wave motion. The barrier was held in place by 'O! ring 
meterial pressed into 0.007 inch slots in the cylinder and 
test module floor. 

The cantilever beams were used to support the cylinder, 
as shown in Figures 6 and 7. Bulkheads were fixed in the 
cylinder at approximately two inches from each end, and the 
fixed ends of the cantilever beams bolted on these bulkheads. 
The free ends of the two beams protruded slightly beyond the 
end of the cylinder and into the plexiglass test module walls 
where they were supported by small self-aligning ball bearings 
pressed into the plexiglass. Both beams were fitted with 
strain gauges and waterproofed using BLH Barrier 'C', as 


shown in Figure 8. 
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Fig. 6 HORIZONTAL FORCE CANTILEVER 
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One beam, positioned with its largest cross sectional 
Gimension vertical, was used to measure Horizontal orcemanad 
the other, with its largest cross section positioned hori- 
zontally, was used to measure vertical force. Calibration 
tests indicated that cross coupling was negligible. 

Dimensions of the beams were chosen to provide sufficient 
flexibility to allow measurement of forces using strain gauges 
and sufficient stiffness to provide a natural frequency that 
was large in comparison with the excitation frequency. In 
addition, the minimum width of the beams was constrained by 
the width necessary to mount the strain gauges. A reasonable 
trade-off was arrived at by making use of the maximum sensi- 
tivity of the amplifier/recorder and determining a maximum 
allowable value for the cantilever beam length. 

To prevent rotation of the cylinder, an offset arm was 
MmeecevQ connect the horizontal force cantilever to the test 
module. A Farber Bearing, AMS 1K7, was used as a wheel at 


the test module end of the arm to reduce friction. 


C. WAVE HEIGHT PROBE 

-A capacitance type probe was chosen to measure wave 
heights in this study. The probe essentially consists of 
2 single insulated wire which acted Po One plate of a capac- 
itor, the water acting as the second plate. The capacitance 
varied as the water level rose and fell. A three-eighths 
inch diameter acrylic rod was attached to the bottom of the 


foil-like support member and the sensing wire was attached 
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at the bottom of the acrylic rod and insulated = theca ae 
tip of the acrylic rod was removable to provide easy changing 
of the sensing wire. The wire was connected at the top to 
a cable connector mounted on a nylon block. Number 30 A.W.G. 
wire with poly thermaleze insulation was used. 

The schematic of the electronic circuit utilized with 
Geis probe is shown in Figure 9. In this circuit, the 
square wave output of a multi-vibrator is used to drive a 
capacitance bridge, one leg of which is varied by the wave 
height probe. The resulting AC error signal is converted to 
abv.C. output signal by the full wave rectifier. The output 
Signal was found to be insensitive to small multivibrator 
frequency changes and linear through the range of the wave 
height probe. A model LD 5211-13 Brush D.C. amplifier was 


Weed With the probe circuit. 


ieee leol PROCEDURE 

Calibration of the probe was conducted with the tank 
filled to the desired level and the test module and wave 
height probe in position. With the probe immersed to a depth 
of nine inches, the amplifier gain was set at desired sensi- 
tivity and ae recorder pen bias adjusted to position the 
recorder to mid-scale. The probe immersion was varied over 
a seven inch range in half inch increments, using a traverse 
mechanism. This initial calibration was checked at severad 


points before and after each set of runs. 
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Fig.9 WAVE HEIGHT PROBE CIRCUIT 
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The force cantilevers were calibrated by loading the 
cylinder with a series of weights, using a pulley arrange- 
ment shown in Figure 10. The load was transmitted to the 
exminder through a series of tapped holes located around the 
Semcumierence at positions horizontally fore and aft, and 
Mertical. The Brush recorder output was calibrated by loading 
Ghe cylinder in half-pound increments. This calibration 
was also checked prior to and after each series of runs. 

After calibration was completed, the wave generator speed 
was set for the desired wave length and a data run commenced. 
During each run, wave generator eccentricity was varied from 
One inch to eight inches, or to a point where the waves 


broke in the channel. 
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FIGURE 10 CANTILEVER CALIBRATION ARRANGEMENT 
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IV. PRESENTATION OF RESULTS AND CONCLUSIONS 


As discussed in the theoretical analysis, dimensionless 
wave force coefficients can be represented as functions of 
the dimensionless water depth (h/a), the dimensionless wave 
height (H/2a), and dimensionless period (gT*/h) parameters. 
feeeou program was developed which, for a given run, held 
water depth and period constant while varying wave height, 
jiaeorder to determine the effect of the dimensionless param-— 
“eeoeon tne horizontal and vertical force coefficients. 
Four series of experimental runs were conducted at relative 
water depths of 9.0, 7.0, 5.5, 4.0. 

The dimensionless wave period parameter was varied for 
each series over a range of 100 to the maximum value attain- 
able with the experimental apparatus. The upper limit in 
wave height was usually provided by breaking. 

A series of representative traces are shown in Figures 
11, 12, 13 and 14 due to the difficulty in presenting all 
the traces. These traces represent experimental data for 
three dimensionless periods from each of the four dimension- 
less water depths. The uppermost trace in each run repre- 
Ssentvs wave height, the middle trace represents vertical 
force and the bottom trace represents horizontal force. The 
experimental eemon the force coefficients for each 
experimental case are given in Appendix B. 

Using the dimensionless period, water depth and waveheight 


as inputs, the computer program was used to generate 
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eorresponding theoretical values of the force vceeti ices 
for each experimental case. This provided a comparison 
between experimental and theoretical values for each run. 
Plots of experimental and theoretical values for maximum 
Hermizontal and vertical force coefficients as a function of 
wave height fOr cacm perted and waver depch svudicdsaree 
included in Appendix C. 

Experimental results and theoretical computations 
indicate that, at the lower values of dimensionless period 
and wave height studied, the wave may be characterized as a 
linear wave. The horizontal force was found to vary sinu- 
soidally, and increased linearly with waveheight. The 
vertical force varied approximately as a negative cosine 
function in this region, having its maximum value in the wave 
trough. The magnitude of the vertical force coefficient 
was found to vary approximately as the waveheight squared. 
ieaewe characteristics are shown most clearly in runs Number 
44 through 47, 71 through 74 and 93 through 95 in Appendix C. 

AS the relative period and wave height increased, the 
wave displayed the sharper peaks and longer, shallower troughs 
that are characteristic of non-linear gravity waves. The 
Iierazontal force, proportional to the fluid particle acceler-— 
ation, tended to reach maximum values closer to the wave 
Peaks, but still remained symmetrical, with respect to the 
wave crest, as indicated by the theoretical traces. The 
theoretical analysis predicted a phase lag between the wave 
crest and the maximum force as small as 24° at the higher 


wave heipnts. 
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The vertical force component began to display positive 
values occurring at twice the frequency of the wave peaks, 
as shown in the experimental traces, Figures 11, 12, 13, 14. 
This is explained by the increasing importance of the lift 
Gemponent of vertical force, which, theoretically, is similar 
to a cosine squared function. At extremely long periods, the 
ieec.torce is clearly dominant, and the maximum force no 
Memeer occurs inethe middle of the wave trough. The»maxiimum 
vertical force occurs when the cylinder is under the wave 
crests where the velocity is greatest. 

At large periods and wave heights, a point is reached 
where flow separation eyes arewno bhenger necisieibive. Thats 
is seen most clearly in the vertical force component, whose 
Velue becomes suddenly lower than the theory predicts, though 
still varying as the waveheight squared. This sudden drop 
in lift is apparently caused by a corresponding decrease in 
the lift coefficient resulting from flow separation. This 
phenomenon is most clearly evident in the graphical presentation 
Gmeruns 51 through 56, and 71 through 73 in Appendix C. 

1). Horizontal Forces 

As previously mentioned, the horizontal forces 
tended to be nearly sinusoidal in nature with a maximum value 
which is quite linear with wave height. It is therefore 
possible to characterize the horizontal force coefficient 
py it's maximum valué (Fy /pea°t). The variation in the 
horizontal force coefficient with dimensionless wave height 


is shown graphically in Appendix C. 
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The experimental values of horizontal force are seen 
to vary linearly with wave height for relative period param- 
eter values from 100 to 120. For relative period values 
greater than 120, the force coefficient varies in an increas- 
ingly non-linear fashion. This is consistent with the theory, 
as shown in Figure 15. However, for the range of period 
parameters studied, the deviation from linearity is relatively 
small. It is therefore possible to display the experimental 
data in terms of the slope of the force coefficient versus 
wave height plots, as shown in Figures 16 and 17. 

The theory predicts the maximum values of horizontal 
force coefficients prior to the wave heights where separation 
effects become apparent. The discrepancy between theoretical 
and experimental values is due to variation of the added mass 
coefficient in a separated flow as well as a net drag force 
mach is no longer zero. 

2). Vertical Forces 

As indicated previously, the vertical force displays 
two regimes; one inertia dominated and one lift dominated. 
ome maximum value of the vertical force coefficient 
(By /pea°t) is seen to vary as the wave height squared 
throughout the range of data taken. This is shown graphically 
in Appendix C. For higher values of relative wave height, 
the magnitude of the vertical force coefficients is less than 
predicted by the theory, due to separation effects, as 


explained above. 
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Ee CONCLUSIONS 

From the experimental and theoretical results, the 
following conclusions are considered warranted: 

in [ne howuZenia ll torees on the cylinder, resi imam 
from the incident gravity waves, increased linearly with 
wave height for relative period values less than 120. For 
relative period values greater than 120, the force coeffi- 
cient was found to vary with wave height in an increasingly 
non-linear fashion. 

2. The vertical forces displayed two regimes; one inertia 
dominated at lower periods and wave heights, and one lift 
dominated at higher values of these parameters. 

3. For high period parameter values, separation effects 
can be expected for relative wave height values above 0.6, 
resulting in a sudden drop in experimental force coefficient 
values from those predicted by the theory. 

4, The theory used in this study provides excellent 
agreement with experimental values for relative period (eT“/h) 
Meawues from 100 to 200, and prior to the onset of flow 


separation. 
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APPENDIX A 


COEFFICIENTS FOR STOKES FIFTH ORDER WAVE THEORY 
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WATER DEPTH IS LISOINGHES.« 


TOTAL NUMBER OF RUNS IS 23. 


Leda weoNy st. 


DIMENSTONLESS PARAMETERS ARE 
WAVE 
HE IGHT 


DATA FOR RUN NUMBER 54. 
WAVE PERIOD IS 
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RUNS IS 23. 


TOTAL NUMBER OF 


DATA FOR RUN NUMBER 55. 


WATER DEPTH IS 11.0INCHES. 


2-647SECONDS. 
DIMENSIONLESS PARAMETERS ARE 


WAVE PERIOD IS 
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TOTAL NUMB=R OF RUNS TS 23. 


DATA FOR RUN NUMBER 56. 


were Cer LS ro Orncres . 


1. 752SECONDS. 
DIMENSIONLESS PARAMETERS ARE 


WAVE PERIOD IS 
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TORAL NUEISER OFURUNS TS "235. 


DATA FOR RUN NUMBER 71. 


WATER DEPTH IS 14.0INCHES. 


2e645SECONDS. 
DIMENSIONLESS PARAMETERS ARE 
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TORE NUMBER OF RUNS IS 23. 


2-365SECONDS. 


DIMENSIONLESS PARAMETERS ARE 


DATA FOR RUN NUMBER720. 


WAVE PERIOD IS 


WATER DEPTH IS 14.0INCHES. 
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TOTAL NUMBER OF RUNS IS 23. 


DATA FOR RUN NUMBER 73. 


WATER DEPTH IS 14.-0INCHES. 


2e170SECONDS. 
DIMENSIONLESS PARBMETERS wane 
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TOUR NO@mwwen OF RUNS TS 25. 


DATA FOR RUN NUMBER 74. 


Wever Berri [S 14 .01N@HES. 


2.022SECONDS. 
NIMENSIONLESS PARAMETERS ARE NON 
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DATA FOR RUN NUMBER 75. 
2e447SECONDS. 
NIMENSTONLESS PARAMETERS ARE 


WAVE PERIOD IS 


TOTAL NUMBER OF RUNS IS 23. 


WATER DEPTH IS 14.0INCHES. 
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OMOOC9ONHNNSSSeSrMNOY 
®©eteteeetet @e © @ © @ @ 


ODODOVODOVIVOOVCNO0O0O0°O 


OFOWTNADAMNOMOMOM 
DAAIN MMOS UV OF Cn O 
ee @e¢ete¢8 #8 &® @ © @ @ 


OVDVODIVVDVDOOO9000O 


OINMFAODODOININOINUDINU 
OMAINM St SLNOM ODORAN 
ee @eee?eee%*#eee @® @ 


ODDVDWOVOOC CO OR Ae 


AO ODS O OOM TOM 
DOVOORAAMNSTTDODOOWNSTWN 
ee ee fe ¢@ @ © @ @ @ @ @ 


DODDONOVDVVV OO FHA 


OOTINISFONNAGWOF- ONO 
eeeeeetetee ®@ @ © © @ @ 


QOOVDVOV OFFA AA ARAN 


DAQ\UMASNURMMODOOOO0O0O 
MODOSFODANMNASTOWNVO 


e*eetrtee © © 0 &© @ © &© @ @ 


DOO HAR ANNOAMOS SN 


rAICASFLVOM ODOMINM SUK 
Oolel@lele@l@lel@ler_ Lal) 
LAVIDY EOLA LAUD LA LALO LA AV LAN) LN ny 
P= Bee RE RE EE EE p- 


13 





TOTAL NOMBER OF RUNS IS 25. 


DATA FOR RUN NUMBER 91. 


WATER DEPTH IS 18.0INCHES. 


2¢665SECONDS. 
DIMENSTONLESS PARAMETERS ARE 


WAVE PERIOD 15 
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TOTAL NUMBER OF RUNS IS 23. 


DATA FOR RUN NUMBER 92. 


WATER DEPTH IS 18.0INCHES. 


2-342SECONDS. 
DIMENSIONLESS PARAMETERS ARE 


WAVE PERIOD IS 
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TOTAL NUMBER OF RUNS IS 23. 


DATA FOR RUN NUMBER 93. 


WATER DEPTH IS 18.0INCHES. 


2-430SECONDS. 
DIMENSIONLESS PARAMETERS ARE 
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TOTAL NUMBER CF RUNS [IS 23. 


DATA FOR RUN NUMBER 94. 


WATER” DEPP IS 1860 INCHES. 


2. OSTSECONGS. 
DIMENSTONLESS PARAMETERS ARE 


WAVE PERIOD IS 
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TOTAL NUMBER OF RUNS [TS 23. 


DATA FOR RUN NUMBER 95. 


WATER DEPTH IS 18.0INCHES. 


1-957 SECONDS: 


DIMENSTONLESS PARAMETERS ARE 
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APPENDIX C. WAVE HEIGHT-FORCE PLOTS 
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APPENDIX D 


COMPUTER PROGRAM 
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